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Abstract 

We study phenomenological implications of a radiative inverse seesaw dark matter model. In 
this model because neutrino masses are generated at two loop level with inverse seesaw, the new 
physics mass scale can be as low as a few hundred GeV and the model also naturally contain dark 
matter candidate. The Yukawa couplings linking the SM leptons and new particles can be large 
leading to large lepton flavour violating effects. We find that future experimental data on \i — )■ cy 
and /x — e conversion can further test the model. The new charged particles can affect significantly 
the h — )■ 77 branching ratio in the SM which is able to explain the deviation between the SM 
prediction and the LHC data. We also study some LHC signatures of the new particles in the 
model. 
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I. INTRODUCTION 

Seesaw mechanism is one of the popular mechanisms[l-3] beyond the standard model 
(SM) which can provide some explanations why neutrino masses are small. Usually the 
seesaw scale is large making LHC study of the new physics scale difficult. The inverse 
seesaw mechanism J] can lower the seesaw scale because in this type of models the light 
neutrino masses are suppressed by higher powers of new scale beyond the SM. If the inverse 
seesaw mechanism is also achieved by radiative correction, the new scale can be even lower. 
Such low new physics scale can lead to large testable effects in various experiments. Recently 
models of this type have been proposed in which inverse seesaw mechanism is radiatively 
realized at two loop level [5] allowing the new physics scale to be in the hundreds GeV range. 
To forbid tree and one loop level neutrino mass generation, new unbroken symmetries are 
introduced. The lightest new particles transforming non-trivially under the new symmetries 
are stable and can play the role of dark matter needed to explain about 23% of the energy 
budget of our universe {(J. 

In this paper we further study some phenomenologies in one of the promising models. 
The model we will study is the U(1)d model discussed in Ref.js|. There are several new 
particles in this model. In this model, the Yukawa couplings linking the SM leptons 
and new particles can be large leading to large lepton flavor violating(LFV) effects. 
We find future experimental data on /x — > and fi — e conversion can further test 
the model. The new charged particles can affect significantly the h — > 77 branching 
ratio in the SM. We find that the new contributions may be able to explain the devia- 
tion between the SM prediction and the LHC data. We provide some details in the following. 



II. THE MODEL 

The model is based on the SU {2) L x U(l)y SM electroweak gauge group with an unbroken 
global U(1)d symmetry. The SM particles do not transform under the U(1)d symmetry. 
New particles in this models are vectorlike leptonic SU {2) L doublets D L R , two scalar singlets 
S, a and a scalar SU(2)l triplet A. Their SM and U(1)d charges are as follows 

D L>R :(2,-1/2)(1), 5:(1,0)(-1), a:(l,0)(2), A:(3,-l)(2). (1) 



2 



In the above the two numbers in the first and the second brackets are the SU(2) L x U(l)y, 
and the U(1)d quantum numbers, respectively. 

The renormalizable terms for Yukawa couplings L D consistent with the symmetries of 
the model are 

L D = -L l YdD r S - D L MD R - l -D L Y L D c L A - ^D C R Y R D R A ] + h.c. (2) 

The allowed renormalizable terms in the potential Vd are given by 

V D = -fJ. 2 H rfH + Xh(H^H) 2 + n 2 s S*S + Xs(S^S) 2 + fcJo + X a (a*a) 2 

+ fJ? A A*A + A^.(A t AA t A) Q + ^ X^ij^j + (pis*S 2 a + X AaH HAa*H + h.c), (3) 

ij 

where the sum £V 1S over & U possible % and j, and % to be one of the H, S, a and A. The 
allowed terms are 

X^H^HtiA)? + X Ha (H^Ha^a) + X HS (H^HS^S) 

+A A s(A t A 1 S t S) + X Aa (A*Aa*a) + X aS (^aS^S) . (4) 

In the above the indices a and j3 indicate different ways of forming singlet. They are given 
by 

(A t AA t A) 1 = A^jAIAm , (A t AA t A) 2 = A^A^A,-, 
(A^AH^H), = A* 3 A l3 H* k H k , (tiAH*H) 2 = A^A^H, (5) 

If both S and A develop non-zero vev's, the Lagrangian will give the usual inverse 
seesaw masses to neutrinos. In that case there will be a Goldstone boson due to breaking of 
the global U(l) d symmetry which may be problematic. To avoid the appearance of massless 
Goldstone boson in the theory, a possible approach is to keep the global symmetry to be 
exact and therefore no Goldstone boson emerges by letting fx 2 to be all larger than zero. This 
also forbids the light neutrinos to have non-zero masses at tree level. However, Majorana 
neutrino masses can be generated at two loop level through the Feynman diagram shown in 
Fig.l. 

Carrying out the loop integrals, one obtains neutrino mass matrix m v in the bases where 
M is diagonal 



/ <H> I <H>\ 

I X- - - f - - X \ 

/ A * \ 

■ 1 ■ X ■ 1 ■ X ■ ' ■ 
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FIG. 1: Two loop Feynman diagram for neutrino mass generation, 
where K k i is defined as 

K k i = ^ 2 (47r) 4 (1 -m 2 /M 2 y ^ m ^ ms ' m ^ ~ 9(m^, M kk ,m s ) 

-gim^ms, M kk ) + g(m^, M kk , M kk )] . (7) 



g(m 1 ,m 2 ,m 3 ) 



f 1 n 2 

/ dx[l + Sp(l - /j 2 ) - ^log/x 2 

Jo 1- fx 2 



with u 2 = ax+b ( 1 x ^ a = ^%,b = Sp(z) is the Spence function or the dilogarithm 

function 

Sp{z) = -p^A dt (8) 

In the above we have assumed that a and the neutral component of A have almost equal 
mass rrifa. 

There are candidates for dark matter in this model. The neutral heavy particles in Dl^r 
and A have non-zero hypercharge and have problems to play the role of dark matter. The 
natural dark matter candidate field is S. It does not have a non-zero hypercharge and does 
not mix with any particle having hypercharge (a mixes with A). As long as dark matter 
properties are concerned, this model is very similar to the real singlet (darkon) model[3] 
and therefore has similar dark matter properties js, 9] and is identical to the complex scalar 



singlet model [lOj with degenerate mass for the real and imaginary parts of S. The term 
important is S^SH^H for dark matter relic density and direct detection. 

The Higgs boson h properties, its mass and its couplings to SM particles (fermions and 
gauge bosons), are the same as those in the SM at the tree level. The recent LHC data 



indicate that the mass is about 126 GeVjll] which can be applied to this model. It has been 
shown that the dark matter relic density and direct detection constraints can be simulta- 
neously satisfied with appropriate dark matter mass. The range of a few tens of GeV for 



dark matter mass is in trouble. However, dark matter mass about half of the Higgs mass or 
larger than 130 GeV is allowed j^j]. In our later discussions, we will take ms = 150 GeV for 
illustration. 

III. NEUTRINO MASSES AND LFV 

The formula in Eq.([7|) determines whether the model is consistent with current data 
on neutrino mixing and masses [l2]. In order to have at least two neutrinos with non-zero 
mass, at least two generations of Dl,r are needed. We will assume that there are three 
of them. The mixing pattern is determined by two Yukawa couplings, Yd and With 
three -Dl,_r, they both are 3x3 matrices. In our numerical calculations, we will assume 
that the flavor structure is determined by the Yukawa coupling Y D with Y D = DdUpmns^d 
with Yi diagonal for both normal and inverted hierarchies for neutrino masses. In our later 
calculations we will use the central values from recent global fit data in Ref. [3] for neutrino 
mixing angles and mass squared differences for both normal and inverted hierarchies (NH 
and IH) for our discussions 

sin 2 fl 12 = 0.307^6(NH,IH);sin 2 # 23 = .386^(NH), 0.392^(IH); 
sin 2 Ota = 0.0241 ± 0.0025(NH), 0.0244±°;°° 23 (IH); 
Sm 2 =m\-m\ = (7.54^;!) x HT 5 eV 2 (NH, IH); 

|Am 2 | = |m 2 - (m 2 + m 2 )/2| = (2.43^f ) x 10" 3 eV 2 (NH), (2.42+°;°D x 10~ 3 eV 2 (IH); 
5 = 194.4°(NH), 196.2°(IH). (9) 

In the following we show two sets of model parameters which can fit known data for 
neutrinos and take them as bench mark values. 

For the normal hierarchy, choosing Y D = diag(l, \A-03, vL77), yo x ^AuH = 10~ 3 , 
Y L = I x 10~ 2 , fji Sa = lOOGeV, m 01 = 300GeV, m s = 150GeV, M u = 500GeV, we can get 
all the three neutrino mass 3.39 x 10~ 2 eV, 3.50 x 10~ 2 eV, 5.98 x 10~ 2 eV, respectively. These 
are consistent with data. 

For inverted hierarchy case, we just need to replace Yd with Yd = 
diag(\/lA6, \/1.48, ^0.100), with all the other parameters unchanged, the neutrino 
masses will be 4.93 x 10~ 2 eV, 5.01 x 10 _2 eV, 3.39 x 10 _3 eV, respectively. Again, these 
numbers are consistent with data. 
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For neutrino masses, the two parameters yo and Xa<tH appear together, but for charged 
lepton LFV processes which happen at one loop level, they only depend on yo- We will 
study how yo is constrained by data from Z, — > ^7 and \i — e conversion. 



FIG. 2: Feynman diagrams for U — > lj"f. 
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FIG. 3: Feynman diagrams for Zj — lj conversion. 

Radiative leptonic decay U — > Ijj can occur at one loop level as shown in FigfSJ By 
attaching 7, and changing 7 into Z, and then let 7 and Z connect to quark, as shown in 
FigJHJ — e conversion can be induced. For our case the Lagrangian responsible for Zj — > Ijj 
and U — lj conversion can be written as 

C = - l 3 o^(A L3i P L + A Rji P R )kF^ + eQ q qYql 3 B L3il ^P L U + H.c] , (10) 

<? 

and the functions A L)R and B L are given by 



A -Y Y* 6 1 F f M *Wi A - mi A 

'S 

Bt,h = YnikYni.^— — — ^ — 77G 



''32TT 2 m 2 s " s rri z s ' J ' rrij 



e 1 „ ,M\ 



2 



„ , . z 2 — 5z — 2 z In z 

Fd(z) = — -, ttt + 



12(2-1) 3 2(Z- l) 4 ' 

^ , . 7^ 3 - 36z 2 + 45^-16 + 6(3^ -2) In z , , 

Gn(z) = : —, . (11) 

v ; 36(1 -z) 4 V ; 

The LFV Zj — > Zj7 decay branching ratio is easily evaluated by 

b & hi) = >4£(i^r + i^H- (12) 

6 



The strength of \i — e conversion is measured by the quantity, B A _+ e = ^convl^capt = ^(a* + 
A(N,Z) -»■ e" + A(N,Z))/T(fi- + A(N,Z) + A(JV + 1,Z - 1)). To obtain the 

conversion rate, one needs to convert the quarks in Eq.flTTl) into relevant nuclei. We will use 
the theoretical values compiled in Ref. 13|. We have|l4| 



where 



and 



B(p-> erf) 



1927i 2 T^ apt 



\D{A)f 



(13) 



(14) 



9lv = ^9lv(u) + 9LV(d)-,9 { LV = 9lv(u) + ^9LV(d), 9LV( q ) = -^eQ q vn il B L . (15) 



The parameters D(A), V^ p,n \A) are nuclei dependent quantities. Several of them are given 
in Ref.Q- 

With the bench mark values for the model parameters fixed, the B(li — > Ijj) and // — e 
conversion rate are all dependent on the coupling constant yn- We now discuss the constraint 
on y D . 

Although // — > cy has not been observed, there are stringent constraint on upper limit 
of 2.4 x 10~ 12 on the branching ratio at the 90% c.l. [ijj]. Experimental sensitivity will 
be improved. We take B(fi — >■ erf) = 1 x 10~ 13 [iiJ as the near future improved MEG 
experimental sensitivity to constrain the parameter yo- There are also bounds for the 
process of r — > fi(e)j. The current 95% c.l. experiment bounds are B{r —> ej) < 3.3 x 10~ 8 , 
B(t ->■ yrf) < 4.4 x 10- 8 [17|. 

There are several measurements of fi — e conversion on various nuclei. The best experi- 
mental bound with the 90% c.l. for Au nuclei is given by B^ e < 7 x 10 — 13 [isl - For Au, the 
relevant parameters determined by method I in Ref. 13j are given by: -D(Au) = 0.189, 
V^(Au) = 0.0974, V^{Au) = 0.146 and fl°_> e (Au)_4= 0.0036(l3|. There are several 



planed new experiments, such as Mu2E [19J ] / C OMET 20] for fi — e conversion using Al. 



The sensitivities are expected to reach 10 16 20j. For Al nuclei, the relevant parameters for 
our calculations aregiven by D(M) = 0.0362, V^(A1) = 0.0161, V^(Al) = 0.0173 and 



R° (A1)= 0.0026[13|. 
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FIG. 4: Figures on the top show constraints on from B(fi — > cy) and \x — e conversion rate 
with current bound (dashed line) and future sensitivity (solid line) for normal(solid line) and 
inverted(dashed line) hierarchy. Figures at the bottom show constraints on yu from r — > ej (left) 
and r — > fij (right). 

The constraints on yjj are shown in FigJH We see that the current upper limits from 
fi — > e7 and fi — e conversion using Au can already constrain yo to be less than 0.2 and 
0.4, respectively. Future fi — e conversion experiments can reach a sensitivity of 0.05 on 
yo- The model will be constrained when new data become available. The constraints from 
r — > /i(e)7 are weaker. 

In the above studies, we have taken some bench mark values to have some ideas about the 
possibility of observing LFV effects. Our studies show that it is possible to have large LFV 
effects observable at near future experiment, in particular for fi — e conversion experiments. 
We, however, should note that from such studies it is not possible to rule out the model 
because the allowed parameters can have large or small LFV effect. For example, for neutrino 
mass generation the scale of neutrino mass depends ony D x XavH, but the leading LFV effects 
we studied XavH does not show up directly. Even assuming all other relevant parameters 
are fixed in the model, by adjusting the size of Aao-h one can have different values of y^ 
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to satisfy constraint on yn from LFV processes. There are some other processes which can 
provide additional test for the model. We find that the correlation of h — > 77 and h — > ^Z 
can be a good indicator. We will discuss this later. 

In the above processes, the new physics scale is set by the masses of particles in the D 
doublet. We now briefly discuss LHC signature of these particles. D°D°, D°D~ and D D + 
can be pair produced through Z, W~ , 7 and Z s-channel exchanges with the cross sections 
of order O(10) fb for mass vtld around 500 GeV. 

For D°D° production, since D° decays into uS, the signature is missing energy which 

would be similar to dark matter signature with emitted photon or gluon from the initial 

quark. The final state is thus a high-pr photon /gluon and missing energy. Detections of 

pi pi 

dark matter pair production processes from CMS 2l| and ATLAS [22j have been performed. 
As no excess from SM predictions observed in both experiments, constraints on dark matter 
mass can then be given accordingly for pair production of dark matter candidates 23( . The 
current data cannot rule out D° of order a few hundred GeV. 

D~ will decay into l~S. The pair production of D°D~ through W~ exchange can be 
searched by pp — > l~ + Efr + X. There will be SM background from W~ — > l~v and 
W~Z/W~ — > l~vvv. Additionally, W + W~ production with leptonic decays will also have 
a possibility to be background when one charged lepton (here / + ) is too soft or too forward. 
To optimize the signal from these backgrounds, one can impose p l T cut on charged lepton. 
With 8 TeV energy at the LHC and 20/6 -1 , it is difficult to cut down the background 
to have enough signal events. We find that it is possible to achieve a discovery level at 
5cr for 14 TeV center of mass frame enery with 300/6 -1 . In Table HJ we show the cross 
section with a selective cut of p l T > 120 GeV. We have chosen the cut for p l T so that the 
signal can be established at 5a level statistically. With a higher cut for p l T , one can have 
a higher significance level, but the event number will be smaller. With this cut of p l T the 
signal is slightly eliminated while the backgrounds are effectively suppressed. Note that, in 
the calculation of signal, we have taken the D doublets with almost degenerate masses and 
used the bench mark values given before. In the narrow width approximation, the cross 
section for I charged lepton in the final state is proportional to X^0 / d^ / d*/ YdYd*) for 
degenerate -Dj. The cross section for pp — )■ l~ + E}? + X, therefore, is almost independent 
of yr). The pair production of D~D + will have the signal plus missing energies. Since 
there are two leptons in the final state, the analysis is more involved. We will not discuss 
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this here. The numbers in Table [U show that with an integrated luminosity of SOOftT 1 for 
centre-of-mass energy at 14 TeV, for both normal and inverted hierarchy cases, signal of 5a 
significance can be achieved using pp — >■ l~ + E]t + X . It is interesting to carry out such a 
search. 





Signal [fb](NH and IH) 


B ackgr ound [fb] 




a(pp^D°Sl-) 


a{W" -> l~p) 


a(W~Z/W- -»■ ri/i/P) 


(t(W + W-) 


14 TeV 


11.1 


14.2 


8.67 x 10 6 


345.3 


1856 


p l T >120 GeV 


9.66 


12.3 


1080 


6.78 


26 



TABLE I: Cross sections of signal with mp = 500GeV and nig = 150 GeV and corresponding 
backgrounds. In both signal and backgrounds, charged lepton of e~ and \x~ are included. 



IV. h — > 77 

There are strong indications from LHC that the Higgs particle has been discovered with 
a mass of 126 GeV whose couplings to gauge bosons are consistent with SM Higgs, but with 
an enhanced h — >■ 77 branching ratio. The experimental value [ll| for this channel is 1.8 ±0.5 
(ATLAS) (1.56±0.43(CMS)) times that predicted by the SM. Recently ATLAS has updated 
their result with[24j 1.8 ± 0.3(stat.)I^^5(sys.)i^^4(theory) times the value predicted by the 
SM. The central value is higher than the SM prediction. If confirmed, new physics beyond 
the SM is required to explain it. In the model we are studying, this can be explained by new 
contribution from charged particles in the triplet scalar A with relatively low mass coupled 
to the usual Higgs boson at loop levels. We now discuss how enhancement can be achieved. 
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26j. 



There have been extensive studies for similar triplet scalar contributions to h — > 77 
Our study is more model inspired, the triplet does not have non-zero vev, and also the A 
does not decay into pure SM particles. The LHC signatures for A particles are different 
than other models. 

In the model we are considering, electroweak symmetry breaking is induced by the non- 
zero vev of Higgs doublet H = (h + , (v + h + H)/\/2) T . The charged h + and the neutral 
fields I are "eaten" by W and Z. The h is the physical Higgs field similar to the one in 
SM. Since this is the only field having a non-zero vev in the theory, at the tree level, the 
Higgs h couplings to gauge bosons are the same as those in the SM. The Yukawa couplings 
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to SM fermions also have the same form as those in the SM. At one loop level, deviations 
start to show up. A particularly interesting one is modification for h — > 77 coupling, due to 
the existence of new charged particles A - ' and their non-zero couplings to h. Note that 
the new particles, do not have strong interactions, the process gg — > h is not affected to the 
lowest order. So the model will not alter the production rate of h predicted by the SM to 
the leading order in agreement with data. 

The couplings of h to A - ' come from A^(A^ AH^ H)i t2 after H develops vev. The 
hAA couplings are given by 

L ~ -[X AH {A + A- + A ++ A— ) + X 2 AH (A ++ A + ^A + A~)]vh . (16) 

Combined with contributions from W and top in the loop, the h — > 77 rate is modified 
by a factor i? 77 = T(h — > 77)c/(i)d 

/T(h 77) S m given by 

p n , v 2 1 r^HA+j^jA v 4(A^ A + A^ A ) 2 

(17) 

where Tj = (m^/4m^), z = t, W 7 , A~ and A . N c is the degree freedom of color and Qt 
is the charge of top quark. Aiijw) and Ay^Tt) come from SM W boson and top quark 
contributions. ^4o( t a) comes from new scalars in the model. They are given by 

A (x) = -x~ 2 [x - f(x)];A l/2 (x) = 2x~ 2 [x + (x - l)/(x)]; 
Ai(x) = -x- 2 [2x 2 + 3x + 3(2x - l)/(x)]; 

Eq. fr20|) tells that new contributions to the ratio i? 77 depend on not only the couplings 
X]f A , but also the masses of the charged scalars. The scalar masses depend on several 
parameters. Neglecting the mixing between a and A , the component fields in A masses 
are given by 

m A0 = /i A + ^A^ A w 2 , 

™l- = Pa + 1 X haV 2 + \>?haV 2 , (19) 
m\— = n\ + -A^ A w 2 + -A 2 ^ 2 , 
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To see how the model can enhance the h — > 77 to be consistent with LHC data, we will 
keep the A mass to be m^o = 300 GeV as used in the discussions on the neutrino masses 
and vary \]f A to obtain the new contributions to R 71 . The results are shown in Fig© and 
Fig© 



We also calculated new contributions to h — > 7Z. We confirm the formalisms in Ref 26]. 



Using these formulas, we obtain the ratio of h — > 7Z in the U{1)d model to that of SM as 
K zj -\ L --rzTi —2 a {z a -,a a -) 

Ji Qn ,r III' / 



SM ""A- 

+ 2gZA ~~ A " ( 2 ^ A + A (z a - , A A - ) } 1 2 (20) 



where Zi = 4mf/m\, Aj = 4m 2 /m| and gz aa = (T A — QaS\v)/ s w c w- Ag M comes from SM 
W boson and top quark contributions and Aq comes from new charged scalars in this model. 
They are given by 



. 2 

Asm — - 

v 



2Q t {T* - 2Q t s 



2 



cot 9 w A l {z w , X w ) + N c ~™™ ~™" WJ Aw(zt, Xt) 

sy/Cw 



A (x,y) = h(x,y), 

Ay 2 (x, y) = h (x, y) - I 2 (x, y) , 

= 4(3- tan 2 6 w )I 2 (x,y) + [(1 + 2s" 1 ) tan 2 6 W - (5 + 2 X - 1 )}I l (x, y), 

where T3 is the third component of isospin of top quark, and Ii, I 2 are given by 

- ^) + l^yf^ - ^ + T^yf^ ~ 9{r% 

h(x,y) = --^—[f(x- 1 )-f(y-% 
2(x-y) 

g(x) = V x' 1 — 1 arcsin \fx. 

In the allowed X l HA space, h —> ^Z will be modified significantly. We show the predicted 
scaling factor R lZ = r(A — >■ jZ) U (r )n /T(h — >■ , yZ) SM in Fig© and Fig© 

To enhance the ratio Ry 7 , negative X]f A are required. With fixed m A o, negative X 2 HA 
implies that m A o > m A - > m A — . From Figj5j we can see that with negative X]f A of order 
(9(1), the ATLAS and CMS results on h — > 77 can be reproduced. If one controls the 
magnitude of X]f A as small as possible from perturbetivity consideration, the optimal values 
for X]f A are around —0.8 and —0.6, respectively. With these values, A~ and A masses are 
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given by 284.5 GeV and 268 GeV. More negative A^ A will make the mass of A smaller 
which may be in conflict with LHC data. We should take |A# A | as small as possible. With 
the same parameters, the predicted value for R lZ are shown in FigfTl 



m A o=300GeV, R yy =1.56±0.43(CMS) m A o=300GeV, R yr =1.8±0.5(ATLAS) 




-2-1 1 2 -2-1 1 2 



FIG. 5: Constraints on A^ A and A^ A with m A o = 300GeV. 



m A o=300GeV, R y7 =1.0 




-2-10 1 2 

^AH 



FIG. 6: Constraints on A^ A and A^ A with i? 77 = 1 for m^o = 300GeV. 




FIG. 7: Scaling factor for h — > Z7 with the same parameters for h — > 77. 
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FIG. 8: Scaling factor for h — > Zj with the same parameters for = 1. 



With X 2 HA > 0, the mass hierarchy for the component fields in the triplet is m A o < 
mi- < m A — • In this case, the new contributions may cancel out if A^ A is kept negative. 
We demonstrate this possibility in FigJHJ , where _R 77 is kept to be 1. For this parameter 
space, the predicted R^z is shown in FigJHJ 

In general there is a correlation between h — > 77 and h — > 7Z, that is, enhancement of 
h — > 77 leads to an enhanced h — > jZ. This fact may be used as a test for this model. 
Should an anti-correlation between h — > 77 and h — > will be confirmed, this model will 
be in trouble. But even if h — > 77 agrees with SM prediction, h — > 7Z can be different as 
can be seen in FigJHJ 

1 2 

Our analysis show that in order to explain the possible enhanced h — > 77, negative A^ A 
of order minus one is needed. If the current data at the LHC will be further confirmed, we 

1 2 

need to check if the required negative A^ A are consistent with other constraints. One of the 
constraints is from the stability of Higgs potential. Here we argue that this is not a problem. 

Potential bounded from below concerns potentials at fields taking large values. Let us 
consider terms involving A#, A A 2 , and Al^ in the case where 5* and a fields are absent and 
carry out an similar analysis as in Ref.[27|. At large values of H and A, the potential is 
given by 



Xhx 2 + (A A + \l V )y 2 + (X 1 HA + X 2 HA 0xy , 



(21) 



where x = WH, y = Tr(A^A), 77 = Tr(A^AA^A) / (Tr(A^A)) 2 , and f = 
AA* H)/(WHTr(A> A)). The ranges for 77 and £ are 1/2 ~ 1 and ~ 1, respectively. 
The positivity conditions of the above require the diagonal elements and the determinant 
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of the following matrix to be positive 

M p = I XH (A " A + Ak0/2 I . (22) 

The parameters need to simultaneously satisfy 

\ H >0, A A + A A >0, A A + ^A A >0, 

A„(A A + A A ) > i(A^ A ) 2 , A„(A A + A A ) > i(A^ A + A^ A ) 2 , (23) 
A# (A A + -A A ) > -(A^ A ) 2 , A^(A A + -A A ) > t(A# a + A^ A ) 2 . 

In the above, we have also taken into consideration of the ranges of rj and £. 

A Higgs mass of 125 GeV, implies Xh = 0.13. Our required X]f A of order minus one 
and the above conditions can be satisfied if one chooses both A A 2 to be positive and satisfy 
A A + \X 2 HA > 1/\ H = 7.5 (with A^ A = -1.0). This condition can be easily satisfie d by 
choosing A A 2 to be about 5 which are well below the unitarity bounds on A A 2 of order 4ir 



27|. 



Our model is more complicated because there are also S and a fields. The term pro- 
portional to A Acr # can be chosen to be small and neglected. The corresponding M p matrix 
becomes a 4 x 4 one. The conditions for potential bounded from below require the diago- 
nal elements, the determinant of the matrix, and all determinants of its sub-matrices to be 
positive. These conditions include the ones discussed above, but have some additional ones. 
For our purpose, we need to fix Xh to be 0.13, and A^ A to be around -1.0 to satisfy the 
positive conditions. Since several new independent parameters \s,cr,Ha,HSASAo-AS come into 
play, one is able to find reasonable parameter spaces to satisfy the conditions. For example, 
with As, CT > 0, Ah ct ,as,Ao-,A5 to be zero. If one requires S to play the role of dark matter, 
Xhs should not be zero [5]. There is a large range for Xhs below 0.03 which can satisfy 
dark matter constraint for dark matter mass around half of Higgs mass and larger than 130 
GeV [5j. The positivity of potential at large values of fields can be satisfied. 

Before closing this section, we make come comments about some effects of the A particle 
at the LHC In the case with A^ A < 0, A is the lightest particle in the A triplet. It is stable 
in the scenario where 2m§ mass is larger than A mass, that is, A — > D~ D~ — > l~Sl~S 
is kinematically forbidden. This is the case for the bench mark values we are using. With a 
mass of order a few hundred GeV, A~' can be produced at the LHC with a cross section of 
order about 10 fb. Although it does not decay into SM particles making the direct detection 
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difficult, being a stable heavy charged particle it does leave tracks in the detector which 
have been searched for at the LHC. The current data from LHC still allow mass of order a 
few Gevba 



|. If it turns out that A mass is large enough, and A — > D~D~ — > l~Sl~S 
becomes kinematically possible, then + EJt+ jets is the signal to search. This has small 
SM background and can be searched at the LHC. 

In the case with A|^ A > 0, A is the lightest particle in the A triplet. It can also be 
copiously produced at the LHC because the mass can be as low as a few hundred GeV. 
Search for this particle is similar to search for dark matter which can annihilate into quarks. 
Some of the processes which can provide information about this particle are single photon 
plus missing energy and mono-jet plus missing energy. ATLAS and CMS experiments at the 
LHC have carried out such studies. At this moment the data are not constraining enough 
to rule out the parameter space we are using [21 



-|23|. But as more data become available, 



the model can be constrained more. 



V. CONCLUSIONS 

We have studied some phenomenological consequences of a two loop radiative inverse 
seesaw with dark matter model with an unbroken global U{l)o symmetry. This model 
allows larger Yukawa couplings and low mass of order a hundred GeV charged new particles 
in the triplet scalar A which couple to the Higgs doublet. The large Yukawa couplings can 
lead to large leptonic flavor changing effects in /i — > and /i — e conversion. The current 
data have already constrained the size of the allowed Yukawa couplings. Future improved 
experiments on /i — e conversion can improve the constraint by several orders of magnitude. 
The existence of low mass charged particle in the triplet A make it possible to enhance the 
h — > 77 to explain the deviation between the LHC data and SM prediction. 
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